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ABSTRACT 

An analysis of the impact of seismic and volcanic activity was carried out at se- 
lected astronomical sites, namely the observatories of El Teide (Tenerife, Canary Is- 
lands), Roque de los Muchachos (La Palma, Canary Islands), Mauna Kea (Hawaii) 
and Paranal (Chile) and the candidate site of Cerro Ventarrones (Chile). Hazard as- 
sociated to volcanic activity is low or negligible at all sites, whereas seismic hazard 
is very high in Chile and Hawaii. The lowest geological hazard in both seismic and 
volcanic activity was found at Roque de los Muchachos observatory, in the island of 
La Palma. 
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1 INTRODUCTION 

Some of the best astrophysical observatories in the world, 
namely the Canarian, Chilean and Hawaiian observatories, 
are located within active geological regions. This is not a 
coincidence, since topography modelled by tectonic and/or 
volcanic activity is a main factor controlling the local at- 
mospheric conditions and hence, the sky transparency that 
defines good astronomical sites. The characterization of as- 
tronomical sites is us ually based on the analysis o f the atmo- 
sphere above them l|Garci'a-Lorenzo et al.l 120071 ). However, 
increasingly larger telescopes (10 and 40 metres classes) de- 
mand stable structures and buildings, and hence geologi- 
cal activity becomes an important par ameter to tak e into 
account in astronomical site ranking l|Sarazin|[2003 ). The 
structures of large telescopes have and will have to with- 
stand the effects associated to seismic and/or volcanic ac- 
tivity, but they also have to minimize the loss of operational 
time, recalling the extreme precision in the alignment of me- 
chanical and optical components. The west coast of Hawaii 
experienced in 2006 October 15, a M=6.7 earthquake fol- 
lowed by a M=6.0 aftershock and many smaller aftershocks. 
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There was not significant structural damage at the telescope 
facilities, however the recovery to full science operability at 
Subaru, Keck I and II and Gemini North telescopes took 
several weeks l|Salmonll200"7l ). 

In this work, we present an analysis of seismic and 
volcanic activity at four world-class observatories, namely 
Roque de Los Muchachos (La Palma, Canary Islands, 
28.75°N, 17.89°E, 2400 m.a.s.L), El Teide (Tenerife, Ca- 
nary Islands, 28.3°N, 16.51°E, 2380 m.a.s.L), Paranal (Chile, 
24.62° S, 70.4°E, 2620 m.a.s.L) and Mauna Kea (Hawaii, 
19.82° N, 155.47° E, 4130 m.a.s.L) and the candidate site of 
Cerro Ventarrones (Chile, 24.35°S, 70.2°E, 2200 m.a.s.L). 
We carried out a comparative analysis of the hazard asso- 
ciated to lava flows, volcanic ashfall, seismicity and ground 
deformation on telescope facilities for a return period of 50 
years, assuming this time span corresponds to the approx- 
imate lifetime f or a large tele scope. Flank collapses of is- 
land volcanoes (|McGuirell200^ . involving mass movements 
of tenths to thousands cubic kilometres of rocks, were not 
considered in this analysis, since their recurrence interval ex- 
ceedes 100000 years for th e Canarian and Hawaiian Islands 
l|Whelan fc Kelletatll2003l ). Moreover, the triggering mecha- 
nisms of these events and the nature of the mass movements 
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are a matter of debate, recalling that humans have not yet 
witnessed any of these catastrophic events. 



2 GEOLOGICAL CONTEXT 
2.1 Tenerife 

Tenerife is the largest island of the Canarian Archipelago 
and one of the largest volcanic islands in the world. It is 
located between latitudes 28-29°N and longitudes 16-17°W, 
280 km distant from the African coast. The morphology of 
Tenerife (see figures 1 and 2) is the result of a complex ge- 
ological evolution: the emerged part of the island was orig- 
inally constructed by fissural eruptions that occurred be- 
tween 12 and 3.3 Ma (|Ancochea et al.lll990l V In the cen- 
tral part of the island, where the observatory is located, the 
emission of basalts and differentiated volcanics gave rise to 
a large central volcanic complex, the Las Canadas Edifice, 
that culminated in the formation of a large elliptical depres- 
sion measuring 16 x9 km^, known as Las Canadas Caldera 
l|Martf et al.l 1199*3 ). In the northern sector of the caldera, 
the Teide-Pico Viejo complex was constructed as the prod- 
uct of the most recent phase of central volcanism. Teide- 
Pico Viejo is a large stratovolcano that has grown during 
the last 175 Ka. The basaltic activity, which overlaps the 
Las Caiiadas Edifice, is mainly found on two ridges (NE 
and NW), which converge on the central part of the island. 
Large scale lateral collapses, involving mass movements of 
hundreds of cubic kilometr es of rock, are respon sible for the 
formation of three valleys (|Masson et al.l [20021 ). namely La 
Orotava (0.35-0.65 Ma), Giiimar (0.75-0.85 Ma) and Icod 
(0.17 Ma) valleys. Eruptive activity at Mount Izana, where 
El Teide observatory is located, ceased more than 300 ka, a f- 
ter the formation of La Orotava valley l|Edgar et al.ll2007l ). 
whereas recorded eruptive activity in Tenerife has consisted 
of six strombolian eruptions, namely Siete Fuentes (1704), 
Fasnia (1705), Arafo (1705), Arenas Negras (1706), Cha- 
horra (1798) and Chinyero (1909). The last three eruptions 
occurred at the NW ridge, the most active area of the is- 
land, together wi th El Teide-Pi co Viejo complex, for the 
last 50,000 years (jCarracedo et al..,2003 . ,2007, ). 



2.2 La Palma 

The island of La Palma is the fifth in extension (706 km^) 
of the Canary Islands and the second in elevation (2426 
m.a.s.l.) after Tenerife (figures 1 and 3). The two main 
stages of the development of oceanic volcanoes, the sub- 
marine and subaerial stages, outcrop in La Palma, since 
the submarine basement or seamount that was built during 
the Pliocene (4 to 2 Ma) were uplifted up to 1.3 km above 
the p resent sea level |staudigcl ct al. 198^; lAncochea et al.l 
1 19941 ). The subaeri al edifice is conformed by a series of over- 
lapping volcanoes (iNayarro fc Coelldll993 : lAncochea et al.l 



to present) in the southern half of the island. The Cumbre 
Vieja ridge is interpreted to be a volcanic rift zone because 
of the prominent north-south alignment of vents, fissures 
and faults. Morphologically, Cumbre Vieja is an extension 
of the Cumbre Nueva ridge and resembles a separate vol- 
cano from the Northern Shield as a result of an emigration 
of the volcanic activity to the south. Historic eruptions on La 
Palma have lasted between 24 and 84 days and were recorded 
in 1470/14 92 (uncertain), 1585, 1646, 1677- 7 8, 1712, 1949 
and 1971 (|Hernandez-Pacheco fc Valid 1 19821 : lRomero-Rui3 
Il99ll ). The Roque de los Muchachos observatory is located 
atop the Taburiente Edifice, at the rim of the Caldera de 
Taburiente, and hence, eruptive activity ceased more than 
0.4 Ma ago. 



2.3 Hawaii 

The island of Hawaii is the youngest island in a chain 
of volcanoes that stretches about 5600 kilometres across 
the northern Pacific Ocean. The island chain results from 
a magma source that originates deep beneath the crust. 
The ocean crust and lithosphere above the magma source, 
within the Pacific tectonic plate, move to the northwest 
with respect to the deep stationary magma source. Over 
a span of about 70 million years, new island volcanoes 
are formed and older volcanoes are carried away from the 
magma source , erode, an d eventually subside beneath sea 
level ( Lockwo od fc Lipma n 1987). Hawaiian volcanoes erupt 
lavas of distinct chemical compositions during four different 
stages in their evolution and growth: a) pre-shield stage, 
representing between 1-2% of the total volume, b) shield 
stage, representing 97-98% of the total volume, c) post-shield 
and post-caldera stage, representing 1% of the total volume 
and d) rejuvenated stage, representing a lot less than 1% of 
the total volume. The island of Hawaii, the largest of the 
entire Hawaiian chain, consists of five main volcanoes: Ki- 
lauea, Mauna Loa, Mauna Kea, Kohala and Hualalai (Fig. 
4) . Kilauea and Mauna Loa volcanoes are in the shield stage 
and erupt frequently, Hualalai and Mauna Kea volcanoes 
are in the post-shield stage and erupt every few hundred to 
few thousand years, and Kohala is dormant, having passed 
through the post-shield stage. Kohala last erupted about 
120 ka, Mauna Kea about 3,6 ka, and Hualalai in 1800-1801 



1 19941 : ICarracedo et al.l 1 19991 : iGuiUou et al.H200ll ): fl) "the 

Northern Shield of the island includes the Garafla shield 
volcano (1.7 to 1.2 Ma), the Taburiente shield volcano (1.1 
to 0.4 Ma) and the Bejenado edifice (0.55 to 0.49 Ma); the 
Caldera de Taburiente, in the center of the old shield, and 
Cumbre Nueva ridge formed by a combination of large land- 
slides and erosion; (2) The Cumbre Vieja Volcano (0.4 Ma 



((EiieraLlITiii). 



2.4 Central Andes, Chile 

The Andes is a mountain chain which extends along the 
western edge of the South American plate. In this region, the 
Nazca Plate thrusts beneath Sou th America at a rate of ap- 
proximately 64 and 79 mm/year (lAngermann et ahll 19991 ) in 
an east-north-east direction, forming the Peru-Chile trench 
in the ocean's floor (Fig. 5). The Andes span a total length of 
about 9000 km extending from Colombia to the Chile triple 
junction around 46° S, which defines the interaction point 
of the South American, Antarctic and Nazca plates. The 
Nazca plate is divided into several segments, whose main 
differ ence are convergence direc tion and dipping angle of the 
slab (|Barazangi fc Isacks|[l97^ ). In general, normal or steep 
dipping slab segments have angles around 30°, whereas in 
shallow segments this angle does not exceed 15°. The angle 
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Figure 1. Simplified shaded-reliof map of the western Canary Islands. The black filled circles indicate the position of El Tcidc Observatory 
(OT) in the island of Tenerife and Roque de los Muchachos observatory (ORM) in the island of La Palma. Filled circles of difi^erent sizes 
and colours indicate the location, magnitude and depth of the earthquakes registered during the period 1973-2008. 



of the segments also define the presence of volcanic activ- 
ity, since in shallow segments there have not been volcanic 
activity at least since the Miocene. 

Both Paranal and Cerro Ventarrones are located in the 
Central Andes region (15°S to 27.5°S), above a steep dip- 
ping segment of the Nazca plate, hence volcanic activity is 
present, although t his takes place more tha n 200 km east- 
wards of the sites (jGregorv- WodzickH I2OO0I ) . In this sense, 
volcanic activity ceased in the area of the observatories more 
than 20 million years ago. The volcanic region of the Cen- 
tral Andes is the largest in the world and it also contains the 
highest active volcanoes, such as LluUaillaco volcano (6723 
m.a.s.L). Unlike the low viscosity-high temperature magmas 
of Hawaii, the high viscosity and low temperature of the 
magmas in the Central andes generate explosive eruptions 
that eject large eruptive columns. 



3 SEISMIC HAZARD ANALYSIS 

Seismic hazard is definied as the probabilistic measure of 
ground shaking associated to the recurrence of earthquakes. 
Seismic hazard maps depict the levels of chosen ground mo- 
tions that likely will be exceeded in specified exposure times. 
Hazard asessment for buildings and other structures com- 
monly specifies 10% chance of exceedance of some ground 
motion parameter for an exposure time of 50 years, corre- 
sponding to a return period of 475 years. The operational 
and survival conditions of telescope structures will depend 
on the ground shaking level (seismic hazard) and the seismic 
design of the structure. Thus, two sites with different lev- 
els of seismic hazard will need different seismic designs for 



the structures to reach the same survival conditions. In this 
sense, larger seism ic hazard means increased building costs 
(jTsang et al.ll2008D . 

Assessing seismic hazard is the first step in evaluat- 
ing seismic risk, obtained by convolving seismic hazard with 
local site effects (topography, geology, amplifications due 
to soil conditions, ...) and vulnerability factors (buildings 
and infrastructures, populations, date and time , ...). The 
basic elements of modern probabilistic seismic hazard are: 
(1) compilation of uniform seismic databases for histori- 
cal (pre- 1900), early instrumental (1900-1964) and instru- 
mental (1964-present) periods; (2) implementation of a seis- 
mic source model to describe the spatial-temporal distri- 
bution of earthquakes, using evidence from earthquake cat- 
alogues, seismo-tectonics, paleoseismology, geomorphology, 
active faults, crustal deformation, among others; (3) eval- 
uation of ground shaking as a function of earthquake size 
and distance; (4) computation of the seismic harzard as the 
probability of ocurrence of ground shaking in a given time 
period, to produce maps of seismic hazard. In order to assign 
a common methodology to infer the probabilistic seismic 
hazard to the different sites, the data from the Global Seis- 
mic Hazard Assessment Program (GSHAP) were analyzed. 
GSHAP (http://www^eismo.ethz.ch/GSHAP/global/) was 
launched in 1992 by the International Litosphere Program 
and implemented in the period 1992-1999 to promote a 
global homogeneous approach to seismic hazard evalua- 
tion (Giardini et al.lll999l) . GSHAP hazard maps depicts 
Peak-Ground-Acceleration (PGA) with 10 % chance of ex- 
ceedance in 50 years, as illustrated in figures 6 to 8 and 
summarized in Table 1. Although all the sites are lo- 
cated within active volcanic/tectonic regions, seismic ac- 
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Figure 2. Simplified shaded-reUef map of the island of Tenerife, indicating the most important geological features (see text for details). 
The light blue and red areas indicate the regions affected by recent (last 10 ka) and historical (last 500 years) volcanic activity, respectively. 
The location of El Teide observatory (OT) is indicated by a black filled circle, whereas filled circles of different sizes and colours indicate 
the location, magnitude and depth of the earthquakes registered during the period 1973-2008. 



tivity strongly depends on the geological setting and thus, 
significant differences in seismicity are expected. This is 
clearly illustrated in the distribution of epicenters recorded 
by the National Earthquake Information Center (NEIC), 
http://neic.usgs.gov for the period 1973 to 2008 (figures 



1 to 5). 

Seismic hazard for the Canary Islands is slightly higher 
than that reported by the GSHAP in the basis of recent 
results. In this sense, recent analysis of historical records 
revealed that earthquakes with intensities of VIII MSK 
took place during volcanic eruptions in Cumbre Vieja, La 
Palma and the er uptions of 1704-17 05 closer to Mount 
Izafia in Tenerife (jRomero-Rui j flQQll ) . Moreover, the fo- 
cal mechanism of the M=5.2 earthquake recorded in 1989 
May 9 (the largest earthquake registered in the Canaries) 
and the aftershock distribution agreeded with the strike of 
a submarine fault para llel to the eastern coast of Tenerife 
l|Gonzalez de Valleio et al. 2006). Hence, the activity of this 
fault and the effect of volcanic activity should be included in 
the hazard calculation, increasing the PGA from 0.15 m/s^ 
to 0.56 m/s^ (0.06g) for eastern Tenerife (including El Teide 
observatory) and 0.5 m/s^ (0.05g) for the rest of the islands. 



In any case, following the convention of the GSHAP, seismic 
hazard at the Canary Islands remains at the lowest level. 

The high rate of earthquake occurrence, coupled with 
the potential for large events, places Hawaii among the ar- 
eas of highest seismic hazard (figures 4 and 6) . Earthquakes 
are relatively common within the Hawaiian Islands, partic- 
ularly on the island of Hawaii, as the result of the crustal 
stresses imparted by the volcanic activity, in particular in 
Kilauea and Mauna Loa volcanoes. Moreover, the weight of 
the island has buckled the earth's crust and isostatic adjuste- 
ments often occur to establish gravitational equilibrium. The 
earthquakes of 2006 October 15 were caused by isostatic re- 
sponse to the stresses generated in the crust and mantle by 
the weight of the island. 

The seismic activity near the Chilean sites of Paranal 
and Cerro Ventarrones (figures 5 and 7) results from the re- 
lease of stresses generated by the subduction of the oceanic 
Nazca plate beneath the South American plate. The fo- 
cal mechanisms of these earthquakes indicate subduction- 
related thrusting, likely on the interface between these two 
plates, as illustrated by the relationship between the loca- 
tion of the hypocentres and the distance from the Peru- 
Chile trench (the longer the distance the deeper the loca- 
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Figure 3. Simplified sliaded-relief map of tlie island of La Palma, indicating the most important geological features (see text for details). 
The light red area indicates the region affected by recent and historical activity of Cumbre Vieja volcano. The location of Roque de los 
Muchachos observatory (ORM) is indicated by a black filled circle. There are not earthquakes with magnitudes M>3 registered in the 
area of La Palma during the period 1973-2008. 




Figure 4. Simplified shaded-relief map of the island of Hawaii, indicating the location of the main volcanoes of the island, namely 
Kilauea, Mauna Loa, Mauna Kea, Kohala and Hualalai. The location of Mauna Kea observatory is indicated by a black filled circle, 
whereas filled circles of different sizes and colours indicate the location, magnitude and depth of the earthquakes registered during the 
period 1973-2008. 
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Figure 5. shaded-relief map of the Central Andes, indicating the most important geological features (see text for details), namely the 
Peru-Chile trench (white solid line) and the volcanic arc (filled black triangles). The location of Paranal observatory and the candidate 
site of Cerro Ventarrones are indicated by black filled circles, whereas filled circles of different sizes and colours indicate the location, 
magnitude and depth of the earthquakes registered during the period 1973-2008. 



Table 1. Seismic hazard expressed in terms of the Peak Ground 
Acceleration (PGA) with 10% chance of exceedance for an expo- 
sure time of 50 years. 



Observatory PGA (g) 



Mauna Kea 0.5 

Paranal 0.47 

Ventarrones 0.42 

Roque de Los Muchachos 0.05 

El Teide 0.06 



tion of the hypocentres) . There is a seismic gap (there have 
not been recorded large earthquakes) between latitudes 18° S 
and 22S°S (see Fig. 5) that mark the zones of the 1868 Au- 
gust 14 (M=9) southern Pe ru and the 187 7 May 9 (M=9) 
northern Chile earthquakes (lKellehedll973 ). Since the esti- 
mated recurrence time for large e arthquakes ranges between 
111 to 264 years (|Nishenkolll99ll '). these regions are likely to 
experience a large earthquake in the future. 



4 VOLCANIC HAZARD ANALYSIS 

In order to assess the volcanic hazard that might face a site, 
it is necessary to know the different types of volcanic activ- 
ity that might be expected and how they present a hazard. 
Volcanic events may be grouped into two types: on one hand 
those directly related with eruptions, like lava flows, pyro- 
clastic flows, pyroclastic airfalls and gas emissions; on the 
other hand those events indirectly related to the eruption 
like landslides, volcanic debris flows (lahars), earthquakes, 
tsunamis and acid rain. To model the efl'ect of near future 
eruptions on astronomical sites, only two volcanic phenom- 
ena have been considered, namely ashfalls and lava flows, 
since they might be considered the most likely volcanic haz- 
ards affecting al l sites, either at s hort (years) or medium 
(decades) terms (|Arafia et al.ll200ol ). 

4.1 Lava flow hazard analysis 

Lava flows are generated by the effusive extrusion of molten 
rock (magma) that pours onto the Earth's surface. The 
hazard from lava flows is the damage or destruction by 
burying, crushing or burning everything in their paths 
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Figure 6. Geological activity map for the island of Hawaii, including Mauna Kea observatory. Green areas indicate the regions where 
there is at least 1 arcsec of ground tilt induced by the fault dislocation represented by the thick solid black line. Areas that might be 
covered by at least 0.01 cm of ashfall after an eruption in the Kilauea volcano, are shown in light orange (see text for details). 



IIMullineaux fc The final morphology of the 

lava flows (length, width, thickness and surace features) is 
determined by various factors such as the rate of effusion, 
the slope of the surface onto which is erupted and viscos- 
ity and chemical composition of the lava. In general, the 
higher the content in silica of the lava the more viscuous it 
becomes. Low silica basaltic lava (<52% Si02) can form fast- 
moving (15-50 km/h) narrow streams or spread out in broad 
sheets up to several kilometres wide. In contrast, higher sil- 
ica lavas (56-70% Si02) tend to be thick, move slowly and 
travel short distances. More silicic lavas (65-75% 8102) tend 
to form mound-shaped features called domes that rarely ex- 
ceed 5 km in diameter. Taking into account that lava flows 
follow valleys and depressions, once potential or actual vents 
are identified, their likely path can be predicted based on lo- 
cal topography. 

In the case of the Chilean sites of Paranal and Cerro 
Ventarrones, the closest volcanoes of the Central Andes are 
at least 200 km distant, being significative longer than the 
typical length of lava fiows in the region (30 km). More- 
over, volcanic activity at the observatory sites ceased more 
than 20 million years ago. Hence, hazard due to lava fiows 
is negligible at Paranal and Cerro Ventarrones. 

The Canarian site of Roque de Los Muchachos is em- 
placed atop the extinct Taburiente volcano, being approx- 
imately 1000 metres above and more than 15 kilometres 
distant from Montana Q uemada, the closest eruptive vent 
of Cumbre Vieja vo lcano (|Hernandez-Pacheco fc Vallslll983 : 
lRomero-Ruiz|[l99ll ). In this sense, hazard due to lava fiows 
is also negligible at Roque de los Muchachos observatory. 



El Teide observatory is emplaced atop Mount Izaiia, where 
volcanic activity ceased more than 300 ka. However, the ob- 
servatory is closer to El Teide-Pico Viejo complex, where 
there has been an intense volcanic activity in the last 150 
ka, whereas several historical effusive basaltic eruptions took 
place at the beginning of the XVIII century within 10 kilo- 
metres from the observatory (Fig. 2). In any case, the caldera 
walls protect the observatory from lava flows originated 
within the caldera area, while local topography and the alti- 
tude of the observatory makes the hazard to lava flows from 
the NE ridge systeni negligible l|Eff-Darwich et all r2008al : 
ICarracedo et aL]|2004l ). 

There is not lava flow hazard at Mauna Kea observa- 
tory due to the eruptive activity of other volcanoes such as 
Mauna Loa, Kilauea or Hualalai volcanoes, as the result of 
the topographic protection given by the altitude the observa- 
tory is located, above 4100 metres above sea level. However, 
Mauna Kea is still an active volcano which last erupted 3600 
years ago and hence, lava flow hazard is low but larger than 
in Chile and the Canaries. 

Following the convention deflned by the United States 
Geological Survey (USGS), the level of hazard for lava fiows 
has been scaled from 1 to 9. being 1 the most hazardous 
(jMuUineaux fc Peters"onlll974l ). The USGS assigns a hazard 
level 7 to the area of Mauna Kea observatory, meaning that 
at least 20 % of this area has been covered by lavas in the 
last 10000 years (see Table 2). Hazard level at El Teide ob- 
servatory was set to 8, since the area covered by lava fiows 
in the last 10000 years is significantly lower than 20 %, Both 
the Chilean sites and Roque de los Muchachos observatory 
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Figure 7. Geological activity map for the Central Andes, including Paranal observatory and the candidate site of Cerro Ventarrones. 
Green areas indicate the regions where there is at least 1 arcsec of ground tilt induced by the fault dislocation represented by the thick 
solid black line (see text for details). Areas that might be covered by at least 0.01 cm of ashfall after VEI=4 and VEI=6 eruptions of 
Lascar volcano are presented in light orange and dark orange, respectively. Peak Ground Accelerations (PGA's) with 10% chance of 
exceedance for an exposure time of 50 years are presented in blue-red contours. 



Table 2. Level of hazard by lava flows that was calculated fol- 
lowing the USGS convention. 



Observatory 


Zone 


Last eruption 


Mauna Kea 


7 


3.6 ka 


Paranal and Ventarrones 


9 


> 1 Ma 


Roque de Los Muchachos 


9 


400 ka 


Teide 


8 


300 ka 



in the island of La Palma have a hazard level of 9, since 
there have not been volcanic activity in the last 60 ka. 

4.2 Volcanic ash risk analysis 

Volcanic eruptions commonly eject into the atmosphere frag- 
ments of molten lava and rocks (pyroclastics and tephra) 
that fall back onto the earth surface. An explosive erup- 
tion can blast this fragments into the air with tremendous 
force. The largest fragments (blocks and bombs) falls back to 
the ground near the vent, usually within 5 km. The smaller 
rock fragments (lapilli and ash) continue rising into the air, 
forming a huge eruption column. Eruption columns can be 
enormous in size and grow rapidly, reaching in some cases 
more than 20 km in height in less than 30 minutes. They 
consist of a lower gas-thrust regions and an upper convec- 
tive region. A column will continue rising convectively until 
its density equals the surrounding atmosphere. When reach- 



ing the troposphere, they will expand laterally and continue 
growing vertically forming a broad umbrella cloud. Once in 
the air, the ash and gases of the column are dispersed by 
the prevailing wind to form an eruption cloud or plume. 
The strength and direction of the wind and the heigth of 
the columns are the main factors controlling the transport 
and final deposition of tephra. 

Ash/tephra fall poses the widest-ranging direct hazard 
from volcanic eruptions. Vast areas (10* to 10^ km^) have 
been covered by more than 10 cm of tephra during some 
large eruptions, whereas fine ash can b e carried out over 
areas of continental size (|Gardeweelll996l ). Burial by tephra 
could collapse roofs, break power and communication lines, 
whereas suspension of fine-grained particles in air affects 
visibility, could damage unprotected machinery, cause short 
circuits in electrical facilities and affect comunications. 

The spatial extend of tephra fall depends on two fac- 
tors, namely the wind distribution and the explosivity of 
the eruption. Wind speeds at different pressure levels for the 
selected sites were coUected from the (NCEP/NCAR) Re- 
analysis database of the National Center for Environmental 
Prediction/National Center for Atmospheric Research. The 
Reanalysis data span from 1980 to 2002, being composed of 
six-hourly and daily U-wind and V-wind components. Wind 
speeds in this database are considered as one of the most 
reliably analysed fields and these data were validated us- 
ing radios onde measurements from the nearest station t o 
each site (|Garcfa-Lorenzo et al.l l2005l : IChueca et al.l |2004| ) . 
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Figure 8. Geological activity map for the western Canary Islands, including El Teide observatory (OT) and Roque de los Muchachos 
observatory (ORM). Green areas indicate the regions where there is at least 1 arcsec of ground tilt induced by the fault dislocation 
represented by the thick solid black line. Areas that might be covered by at least 0.01 cm of ashfall after an eruption in the central 
volcanic complex of Tenerife and Cumbre Vieja volcano in La Palma, are shown in light orange. 



The vertical distribution of wind direction at all sites (fig- 
ures 9 to 11) shows a similar behaviour above the 700 mBar 
pressure level (approximately 2500 m.a.s.L), being the pref- 
erential wind direction eastwards. Hence, wind is blowing 
approximately in the same direction above all the selected 
sites. 

The second major factor affecting the extend of ashfall 
is the explosivity of the eruption, expressed as the Volcanic 
Explosivity Index (VEI). In this sense, a VEI=3 eruption 
produces eruptive columns between 3 to 15 km high above 
the vent, a VEI=4 eruption produces eruptive columns be- 
tween 10 to 25 km high, whereas eruptions with VEI larger 
than 4 are associated to eruptive columns higher than 25 
km. The Explosivity Index for eruptions in the Central An- 
des that have been measured in historical times ranges from 
the VEI =4 su b-plinian eruption of Lascar volcano in 1967 
or 1993 l|Gardeweg 199Q ) to the VEI=6 eruption of Huay- 
naputina volcano in February 1600, corresponding to the 
largest volcanic explosion in South America in historic times. 

The Explosivity Index for eruptions in the island of 
La Palma is low (VEI<3), although eruptive columns up 
to 4 km high have been recorded in historical times when 
the uprising magma encount ered groundwater, p roducing 
phreatomagmatic explosions (| Romero- Ruij Il99lh . In the 
case of the island of Tenerife, explosive eruptions might 
occur in the central part of the island, associated to El 
Teide-Pico Viejo stratovolcano and its peripherical vents. 
Last explosive event took place 2 ka and corresponded to 
the VEI=3 sub-phnian eruption of Montafia Blanca volcano 



dAblav et al.ll 19951 ') . Due to the complex nature of the central 
volcanism in Tenerife, other types of explosive eruption have 
occurred and might occur in the future, namely phre atomag- 
matic and/or violent strombolian eruptions ( Mar tf et al.l 
l2008h . The Explosivity Index for Hawaiian eruptions is simi- 
lar to that found in the Canarian island of La Palma, namely 
VEI<3, but with chances of phreatomagmatic eruptions like 
those reported for the Kilauea volcano 2700, 2000, 1100, 600 
years ago and the er utptions occurred in 1790 and May 1924 
iMastin et al.ll 19981 ). 

Unlike the hazard analysis for lava flows, there is 
not established convention to define the hazard due to 
ash/tephra fall. H ence, a numerical simulator of volcanic 
asfaU, TEPHRA2 jConnor et aLlbOOll ). was used to analyze 
the extend of ash deposits at each site, considering the typi- 
cal wind conditions (figures 9 to 11) and the VEI of the vol- 
canoes located in the surroundings of the sites. In this sense, 
we supposed two different eruptive events (VEI=4 and 6) for 
the Lascar volcano in Chile, a 8 km high eruptive column 
for a VEI=3 event in Montana Blanca (Tenerife) and 4 km 
high basaltic eruptive columns for VEI=2 eruptions in Cum- 
bre Vieja (La Palma) and Kilauea (Hawaii). The results of 
TEPHRA2 (figures 6 to 8) show that El Teide observatory 
could be affected by a significant deposition of ash (up to 
1 cm) due to the prevaling wind conditions, recalling that 
approximately 60 % of the time the wind blows eastwards. 
The Chilean observatory of Paranal and the candidate site 
of Cerro Ventarrones could be affected by a VEI=6 eruption 
if the wind blowed westwards, a condition that happens only 
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Figure 9. Frequency distribution of wind direction (rod lines) 
and velocity (black lines) on the Canary Islands at five different 
pressure levels, namely 925, 700, 400, 200 and 100 mbar. 



Figure 11. As in Fig. 9, but for the case of Paranal and Cerro 
Ventarones, Chile. 
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Figure 10. As in Fig. 9, but for the case of the island of Hawaii. 



0.5 % of the time. In the case of La Palma and Hawaii, the 
chances to have significant depositions of ash at the obser- 
vatory site are negligible, as a result of the distance of the 
sites relative to the eruptive vents and the direction of the 
prevailing wind. 

In an attempt to assign a probabilistic measure of the 
hazard associated to ash/tephra fall, we used t he method- 
ology implemented bv lNewhall fc HoblittI (|2002l ). This tech- 
nique assigns the probability that an area might be cov- 
ered by at least 10 cm of ash during the next 50 years, in 
terms of the VEI and distance of the volcanic eruption, un- 
der the supposition of two different return periods for explo- 
sive events, 200 and 2000 years (Table 3). The probabilities 
of being covered by ash in the astronomical sites of Roque 
de los Muchachos (La Palma) and Mauna Kea (Hawaii) are 
negligible, whereas the probabilities in Chile and Tenerife 
depends on the return period of explosive events. Several 
explosive eruptions (VEI > 3) have occurred in the Central 
Andes in the last 100 years l|Gardeweg||l996l ). whereas only 
two of thes e events have been reported in Tenerife in the 
last 35 ka l|Carracedo et all 120071 : iMarti et al.] 120081 ). The 
low rate of occurrence of explosive eruptions in Tenerife is 
also supported by the fact that there are not tephra deposits 
from the Teide-Pico Viejo complex at the observatory, al- 
though this is located in an area with probabilities of being 
covered by ashfall. 
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Table 3. Probabilities that an area might be covered by at least 
10 cm of ash for an evaluation period of 50 years and two different 
return periods of eruptive events, namely 200 and 2000 years. 



Observatory 


A 


B 


Mauna Kea 


0% 


0% 


Paranal and Ventarrones 


7.5% 


0.75% 


Roque de Los Muchachos 


0% 


0% 


El Teide 


12% 


1.2% 



5 GROUND DEFORMATION 

Present-day telescopes reach precisions in pointing and 
tracking below 1 arcsec, hence stability is required in both 
the telescope structure and the ground. This is partic- 
ularly important for long baseline interferometers, where 
precisions better than 20 /im are require d in the align - 
ment of mechanical and o ptical components (|Scholleill2008h . 
lEff-Darwich et al.l (l2008bl ) studied the effect of ground tilt at 
astronomical sites induced by volcanic reactivation of mag- 
matic systems associated to stratovolcanoes and to dike in- 
trusions. These authors concluded that significant ground 
tilts of at least 1 arcsec are found within 10 to 15 kilome- 
tres from the magmatic system. According to the VEI of 
the eruptions and the distance to the different astronomi- 
cal sites. El Teide observatory could be affected by ground 
tilt larger than 1 arcsec in the case of reactivation of the 
central volcanic system (El Teide-Pico Viejo stratovolcano 
and peripherical edifices) and/or a close- by dike intrusion 
similar to those responsible for the historical eruptions in 
the island. Moreover, Mauna Kea observatory might also be 
affected by significant ground tilt as the result of the re- 
activation of this volcano; however, there is not significant 
ground tilt induced by activity from Kilauea and/or Mauna 
Loa volcanoes. Ground tilts at Roque de los Muchachos ob- 
servatory and the Chilean sites are negligible as the result 
of the distance to the lik ely active vents, as observed by 
iPritchard fc Simons! (|2002l ) in the Central Andes volcanoes. 

We carried out a theoretical analysis to infer the ground 
deformation, in terms of ground til t, associated to a dislo - 
cation induced by a fault. Following iFeieel fc Dupri l|l999h . 
we specified nine parameters describing the rectangular fault 
patch: three centroid coordinates; strike and dip; along- 
strike length L and down-dip width W; and the slip vector 
of the left-lateral Ul and up-dip U2 components. These pa- 
rameters were calculated from the linear relations between 
ear thquake magnitud e M a nd fault morphologies obtained 
by iPapazachos et al] l|2004l ') for different types of tectonic 
settings. Following the probabilistic hazard maps calculated 
by GSHAP, likely earthquakes and the corresponding fault 
parameters were defined at each site (see Table 4). Abso- 
lute vertical displacements and ground tilts associated to 
the fault dislocations were calculated t hrough the RNGCIfN 
software developed bv Feigel fc Dupre (1999). It is observed 
(figures 6 to 8) that significant ground tilt of at least 1 arcsec 
might affect the Chilean sites, whereas the effect of ground 
tilt is less likely in Hawaii and Tenerife and negligible in La 
Palma. 



Table 4. Fault parameters used to calculate ground tilt at the 
different astronomical sites, namely earthquake magnitude M, 
Depth of the top of the fault (in km), dip angle of the fault (de- 
grees), length and width of the fault (in km) and dip displacement, 
U2 (in cm). 



Observatory 


M 


Depth 


Dip 


Length 


Width 


U2 


Mauna Kea 


7.2 


30 


25 


55 


20 


240 


Central Chile 


8.0 


40 


25 


160 


70 


220 


Canary Islands 


6.8 


25 


60 


48 


10 


100 



6 CONCLUSIONS 

An analysis of geological hazards associated to seismic and 
volcanic activity (lava flows and ashfall) was carried out at 
selected astronomical sites, namely El Teide and Roque de 
los Muchachos in the Canary Islands, Mauna Kea in Hawaii 
and the Chilean observatory of Paranal and the candidate 
site of Cerro Ventarrones. A common methodology was used 
to characterize the geological hazard, expressed in terms of 
probabilities of occurrence in the next 50 years, recalling 
that this period of time corresponds to the expected life- 
time of a telescope. Large telescopes have to withstand the 
effects associated to seismic activity, but they also have to 
minimize the loss of operational time, recalling the extreme 
precision in the alignment of mechanical and optical compo- 
nents. The hazard from lava flows is the damage or destruc- 
tion by burying, crushing or burning astronomical facilities 
in their paths. Burial by ash/tephra could collapse roofs, 
break power and communication lines, whereas suspension 
of fine-grained particles in air affects visibility, could dam- 
age unprotected machinery, cause short circuits in electrial 
facilities and affect comunications. Tectonic and/or volcanic 
activity might also induce subtle ground tilt that could also 
result in the disalignment of mechanical and optical compo- 
nents, in particular in the case of interferometric observa- 
tions. 

The lowest geological hazard in both seismic and vol- 
canic activity was found at Roque de los Muchachos obser- 
vatory, in the island of La Palma. Seismic hazard is also 
low at the other Canarian site. El Teide observatory, since 
seismic activity in the Canary Islands is low in both number 
and magnitude of earthquakes. On the contrary, seismic haz- 
ard is very high in Paranal and Ventarrones (Chile) and in 
Mauna Kea (Hawaii). Hazard associated to lava flows dur- 
ing a volcanic eruption is not significant at any site, as the 
result of low volcanic activity in the regions where the sites 
are emplaced, topographical protection or distance to the 
eruptive vents. Hazard associated to volcanic ashfall is neg- 
ligible at Mauna Kea and Roque de los Muchachos and low 
to moderate in Tenerife and the Central Andes, depending 
on the prevailing winds and the still poorly known explosive 
volcanic activity in these regions. Ground tilt induced by 
volcanic activity might be significant at Mauna Kea and El 
Teide observatories, whereas tilt produced by fault disloca- 
tion is significant at the Chilean observatories. 



12 Eff-Darwich et al. 



ACKNOWLEDGMENTS 

This work has made use of the NCEP Reanalysis data 
provided by the National Oceanic and Atmospheric Ad- 
ministration / Cooperative Institute for Research in En- 
vironmental Sciences (NOAA / CIRES) Climate Diagnos- 
tics Center, Boulder, Colorado, USA, from their web site at 
|http://www.cdc.noaa.gov/| 



REFERENCES 

Ablay, G.J., Ernst, G.G.J. , Marti, J., Sparks, R.S.J., 1995, 

Bull. VolcanoL, 57, 337 
Ancochea, E., Fuster, J.M., Ibarrola, E., Cendrero, A., 

Coello, J., Hernan, F., Cantagrel, J.M., Jamond, A., 1990, 

J. VolcanoL Geotherm. Res., 44, 231 
Ancochea, E., Hernan, F., Cendrero, A., Cantagrel, J.M., 

Fuster, J.M., Ibarrola, E., Coello, J., 1994, J. Volcano!. 

Geotherm. Res., 60, 243 
Angermann, D., Klotz, J., Reigber, C, 1999, Earth and 

Planetary Science Letters, 171, 329 
Arafia, V.A., Felpeto, M. Astiz, A. Garcia, A. Ortiz, R., 

Abella, R., 2000, J. VolcanoL Geotherm. Res., 103, 377 
Barazangi, M., Isacks, B., 1976, Geology, 4, 686 
Carracedo, J.C., 1994, J. VolcanoL Geotherm. Res., 60, 225 
Carracedo, J.C., Day, S., Guillou, H., Gravestock, P., 1999, 

GSA Bulletin, 111, 755 
Carracedo, J.C., Paterne, M., Guillou, H., Perez Torrado, 

F.J., Paris, R., Rodriguez Badiola, E., Hansen, A., 2003, 

Estudios Geol. 59, 15 
Carracedo, J.C., Guillou, H., Paterne, M., Scaillet, S., 

Rodriguez Badiola, E., Paris, R., Perez Torrado, F.J., 

Hansen, A., 2004, Estudios Geol. 60, 63 
Carracedo, J.C., Rodriguez-Badiola, E., Guillou, H., Pa- 
terne, M., Scaillet, S., Perez- Torrado, F.J., Paris, R., Fra- 

Paleo, U., Hansen, A., 2007, GSA Bulletin, 119, 1027 
Chueca, S., Garcfa-Lorenzo, B., Mufioz-Tufion, C, Fuen- 

salida, J.J. 2004, MNRAS, 349, 627 
Connor, C.B., B.E. Hill, B. Winfrey, N.M. Franklin, and 

P.C. LaFemina, 2001, Natural Hazards Review, 2, 33 
Edgar, C.J., Wolff, J.A, Olin, P.H., Nichols, H.J., Pittari, 

A., Cas, R.A.F., Reiners, P.W., SpeU, T.L., MartiJ., 2007, 

J. VolcanoL Geotherm. Res., 160, 5985 
Eff-Darwich, A., Garcfa-Lorenzo, B., Bonatto, L., 

Hernandez-Gutierrez, L. E., Vifias, R., Rodriguez- 

Losada, J. A., Blanco, M. J., & Mufioz-Tufion, C, 2008, 

Proceedings of the SPIE, 7012 
Eff-Darwich, A., Grassin, O., Fernandez, J., 2008, Pure and 

Applied Geophysics, 165, 1049 
Feigel, K.L., Dupre, E., 1999, Comp. Geosci, 25, 695 
Frey, F.A., Kennedy, A., Garcia, M.O., West, H.B., Wise, 

W.S., Kwon, S.T., 1984, Eos, Transations, American Geo- 
physical Union 65, 300 
Garcfa-Lorenzo, B., Fuensalida, J. J., Mufioz-Tunon, C, 

Mendizabal, E., 2005, MNRAS, 356, 849 
Garcfa-Lorenzo, B., Fuensalida, J. J., Rodriguez- 
Hernandez, M. A. C, 2007, Proceedings of the SPIE, 

6747 

Gardeweg, M., 2000, MM A Site East of San Pedro De At- 
acama North Chile Volcanic Hazards Assessment and Ge- 
ologic Setting, MMA Memo 251 



Giardini, D., Griinthal, G., Shedlock, K.M., Zhang, P., 

1999, Annah de Geofisica, 42, 1225 
Gonzalez de Vallejo L., Garcfa-Mayordomo J., Insua J.M., 

2006, Bull. Seism. Soc. Am. 84, 2040 
Gregory-Wodzicki, M., 2000, Geol. Soc. Am. Bull., 112, 

1091 

Guillou, H., Carracedo, J.C., Duncan, R., 2001, J. Volcano!. 

Geotherm. Res., 106, 175 
Hernandez-Pacheco, A., Vails, M.C., 1982, Arquipelago, 

111, 83 

Kelleher, J. A., 1972, J. VolcanoL Geotherm. Res., 77, 2087 
Lockwood, J. P., Lipman, P.W., 1987, Volcanism in Hawaii: 

U.S. Geological Survey Professional Paper 1350, 509 
Martf, J., Mitjavila, J., Arafia, V., 1994, Geol. Mag., 131, 

715 

Martf, J., Geyer, A., Andujar, J., Teixido, F., Costa, F., 
2008, J. VolcanoL Geotherm. Res., 178, 529 

Masson, D.G.; Watts, A.B., Gee, M.J.R., Urgeles, R., 
Mitchell, N.C., Le Bas, T.P., Canals, M., 2002, Earth Sci. 
Rev., 57, 1 

Mastin, L., Christiansen, R., Swanson, D., Stauffer, P., 
Hendley II, J., 1998, U.S. Geological Survey Fact Sheet 
132, 98 

MuUineaux, D., Peterson, D., 1974, U.S. Geological Survey 

Open File Report, 72-239 
Navarro, J. M., Coello, J., 1993, Mapa geologico del Parque 

Nacional de la Caldera de Taburiente. ICONA. Ministerio 

de Agricultura, Pesca y Alimentacion 
Nishenko, S.P., 1991, Pure Appl. Geophys., 135, 169 
Newhall, C, Hoblitt, R., 2001, Bull. Volcano!., 64, 320 
Papazachos, B.C., Scordilis, E.M., Panagiotopoulos, D.G., 

Papazacos, C.B., Karakaisis, G.F., 2004, Bulletin of the 

Geological Society of Greece, 36, 1482 
Pritchard, M.E., Simons, M., 2002, Nature, 418, 167 
Romero-Ruiz, C, 1991, Las manifestaciones volcanicas 

historicas del Archipielago Canario, Consejerfa de Polftica 

Territorial, Gobierno de Canarias 



Salmon, D., 2007, http://www.gemini.edu/earthquake/agenda.html 



Sarazin, M., 2002, ASPC, 266, 342 

Scholler, M., 2008, FAS Publications Series, 28, 1 

Staudigel, H., Feraud, G., Giannerini, G., 1986, J. Volcano!. 

Geotherm. Res., 27, 299 
Tsang, D., Austin, G., Gedig, M., Lagally, G., Szeto, K., 

Sagals, G., Stepp, L., 2008, Proceedings of the SPIE, 7012, 

154 

Whelan, F., Kelletat, D., 2003, Progress in Physical Geog- 
raphy, 27, 198 
McGuire, W.J., 2006, Phil. Trans. R. Soc. A, 364, 1889 



